Phase Structure of Repulsive Hard-Core Bosons in a Stacked Triangular Lattice 
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In this paper, we study phase structure of a system of hard-core bosons with a nearest-neighbor 
(NN) repulsive interaction in a stacked triangular lattice. Hamiltonian of the system contains two 
parameters one of which is the hopping amplitude t between NN sites and the other is the NN 
repulsion V . We investigate the system by means of the Monte- Carlo simulations and clarify the 
low and high-temperature phase diagrams. There exist solid states with density of boson p = | and 
|, superfluid, supersolid and phase-separated state. The result is compared with the phase diagram 
of the two-dimensional system in a triangular lattice at vanishing temperature. 
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In recent years, cold-atomic systems are one of the 
most intensively studied topics not only in atomic physics 
but also in condensed matter physics. In particular, 
study on cold-atom systems put on an optical lattice may 
give very important insight about dynamics of strongly- 
correlated particle systems^. For systems in an opti- 
cal lattice, interactions between atoms, dimensionality 
of system, etc. are highly controllable, and effects of im- 
purities are strongly suppressed. Among many interest- 
ing topics, possibility of the appearance of the supersolid 
(SS) has been investigated intensively. In the SS, the 
superfluidity and a diagonal long-range order (DLRO) 
such as density wave coexist. It has been clarified that 
the SS cannot exist in a hard-core (HC) boson system in 
a square lattice unless long-range interactions between 
bosons are included^. On the other hand on a triangular 
lattice, it was shown that the SS can appear as a re- 
sult of the competition between the particle hopping and 
nearest-neighbor (NN) repulsion^i^. 

In this paper, we shall pursue the above problem, 
i.e., we shall consider a three-dimensional (3D) hard-core 
(HC) boson system in a stacked triangular lattice and 
study a finite-temperature (T) phase diagram. As we 
consider the 3D system, there exist finite-temperature 
(T) phase transitions in addition to "quantum phase 
transition" . Hamiltonian of the system is then given by 



+V ^ riiUj - ji'^rii, 

{{ij)} i 



(1) 



where i and j denote site of 3D stacked triangular lat- 
tice, is annihilation operator of the HC boson, Ui is 
the number operator rii = ^j^^ and fi is the chemical 
potential for the grand-canonical ensemble, (z, j) denotes 
the NN sites in the 3D stacked triangular lattice, whereas 
((i, j)) those of the 2D triangular lattice, t is the hopping 
amplitude in the 3D lattice, whereas V{> 0) denotes the 
NN repulsion in the 2D triangular lattice. The HC bo- 
son operators satisfy the following (anti) commutation 



relations. 



[-^i, <Pj] = [4, 4] = [4, <Pj] = 0, for i ^ j, 

{4>i,4} = 1. {4iAi} = 0, {4^4} = 0- 



(2) 



The above HC boson can be expressed in terms of the 
Schwinger boson Wia- {cr = 1,2), 



(3) 



where the physical-state condition of the Schwinger is 
given as. 



(7=1,2 



(4) 



where |Phys)^(; is the physical Hilbert space of Wi^. 

Partition function of the system at temperature (T) 
Z = Trexp(— /3i^) with P = is given as follows, 

Z = j\[[dwi{T)]e^Y>{fdTA{T)), 

^(^) = - E *^^'(^) - ^("^(^))' (5) 

where [dwiir)] = S{J2a '^icrWia - l)cPwi{T). In the 
present study, we shall ignore the r-dependence of Wia-, 
and then focus on the finite-T physical properties, 



Z' = j llidwijexpi-pHiw)). 



(6) 



In the previous papers^*^, we discussed that state with a 
LRO found in the above approximation at finite T sur- 
vives to T = unless it is replaced by another ordered 
state. This expectation will be again confirmed by the 
study in this paper through the comparison between the 
obtained results and those of system in a 2D triangular 
lattice at T = 0^. From Eq.(|6]), we introduce dimen- 
sionless parameters ci = and C2 = PV. To obtain 
phase diagram, we fix the parameter ci and investigate 
the system by varying values of C2 and /i (the grand- 
canonical ensemble), or the density of (the canoni- 
cal ensemble). These two calculations, i.e., the grand- 
canonical and canonical ensembles, are complementary 
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with each other. For the simulations, we employ the 
standard Monte-Carlo Metropolis algorithm with local 
update calculation. The typical sweeps for measurement 
is (30000 - 50000) x (10 samples). 
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FIG. 1. Phase diagram in the grand-canonical ensemble 
for C2 = V/kuT = 50.0, which corresponds to the low- 
temperature phase. 

We first show the phase structure of the system in 
the grand-canonical ensemble. To investigate the phase 
structure, we measured the internal energy the spe- 
cific heat C, the particle density p, etc, which are de- 
fined as = {H)/L^, C={{H- Ef)/L^ and p = {m), 
where L is the linear system size. We numerically studied 
the system at various T's and found that the phase dia- 
gram changes continuously. For the low-T phase, we set 
C2 = 50 and obtain the phase diagram in the t/V — ji/V 
plane. See Fig(TJ There are full and empty states in 
the phase diagram, in which all sites of the system are 
occupied or empty, respectively. For t/V <C 1, there ap- 
pear "solid" states with the density p = \ and p = |. 
It is easily seen that these states in which one of three 
sites is filled (empty) in a a/3 x a/3 ordering have min- 
imum energy for t/y <C 1^. As t/V is increased, there 
appear the SS and superfiuid (SF) states. The obtained 
phase diagram in Fig(T]is essentially the same with that 
of the system in the 2D triangular lattice at T = 0^. The 
three dimensionality of the present system preserves the 
ordered states at T = in 2D up to the low but finite T. 

In the following, we shall exhibit the MC calculations 
that derived the above phase diagram in Figdl First in 
FiglSl we show C and p as a function of p for ci = 5 
and C2 =50, i.e., t/V = 0.1. It is obvious that there is a 
particle-hole symmetry. Singular behavior of C indicates 
that there are totally six phase transitions and the step- 
function like behavior of p indicate that phase transitions 
at p/V 0, 3 and 6 are of first order. 

In order to verify the physical meaning of each phase, 
we measured the particle-density correlation function 
n(r) and the HC boson correlation function G(r), which 
are defined as follows, 



n(r) = G(r) = (^^i+r). 

It is expected that in the solid states with p = 



(7) 
and 



1 

3 

|, n(r) exhibits a specific correlation with 3a periodicity 



FIG. 2. Internal energy, specific heat and particle density as 
a function of p/V for ci = 5 and C2 = 50 {t/V = 0.1) in the 
grand-canonical ensemble. System size L — 24. 
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FIG. 3. Correlation functions n{r) and G{r) in various phases 
for ci = 5 and C2 = 50 {t/V — 0.1) in the grand-canonical 
ensemble. L = 24. 



(a = lattice spacing), whereas vanishing superfiuid cor- 
relation, i.e., G{r) ^ for r ^ large. On the other hand 
in the superfiuid state, n(r) exhibits a homogeneous dis- 
tribution of bosons and there exists a finite superfiuid 
correlation. Supersolid has finite correlation of the both 
n(r) and G(r). We have verified all of the above expec- 
tations by the MC simulations, and here we show some 
of the results in Figl3l In the solid states, the HC bosons 
are localized in a specific pattern to avoid repulsions and 
as a result a density wave appears. Then, the superfiuid 
correlation decays very rapidly there. The superfiuid is 
a homogeneous state with nonvanishing HC boson LRO. 
On the other hand in the SS, both the density wave and 
HC boson LRO appear. 
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FIG. 4. Snapshot of density n(r) in various phases for ci = 5 
and C2 = 50 (tjV — 0.1) in the grand-canonical ensemble. 

p — \ solid, Upper right: p ■ 
SF, and Lower right: SS. 



Upper left: p = | solid, Upper right: p = ^ solid. Lower left: 




FIG. 5. Phase diagram in the grand- canonical ensemble for 
C2 = V/kBT = 10.0 

Snapshots are useful to get an intuitive physical picture 
of each phase. In FiglH we show snapshots of particle 
density of each phase. In the solid states, particles are 
localized, whereas in the superfiuid a homogeneous state 
is realized. In the SS, the high-density region fluctuates 
and localized domains come to overlap with each other. 

It is interesting to see how the phase diagram changes 
as T is raised. As we are studying the 3D system, we 
can address this interesting problem. We show the phase 
diagram of the system with C2 = PV = 10.0 in FigEl 
It is obvious that the parameter regions of the SS and 
also SF are getting small compared to those at low T. 
Phase boundary of the two solids with p = | and | be- 
comes obscure. The calculated specific heat C and the 
particle density p are show in FigE] as a function of the 
chemical potential. For ci = 1.0 and C2 = 10.0, there 
are two sharp peaks in C at p/V 1.0 and 5.5, which 
are related with each other by the particle-hole symme- 
try, and correspond to the empty-solid/solid- full phase 
transitions. There also exists a small V-shape hollow at 
/i/V 3.0, i.e., p 0.5. Calculation of p exhibits rather 
smooth behavior without the sharp step- function like dis- 
continuities that exist at low T. The obtained correlation 
functions in Fig |7] indicate that a mixed state of the two 
solids with p = | and p = | appears for 1.0 < p/V < 5.5 
and the density of each solid depends on value of p/V. 

Let us turn to the phase diagram of the system in the 
canonical ensemble. We use the MC simulations with lo- 
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FIG. 6. Specific heat and density as a function of p/V for 
ci = 1.0 and C2 = 10 (i.e., t/V = 0.1). 
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FIG. 7. Correlations of density and particle for ci 
C2 = 10. 
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cal updates that conserve the total boson number. We 
consider the low-T case with C2 = 50 and show the phase 
diagram in the t/V — p plane in FiglHl We show only 
the region < p < 0.5 because the particle- hole sym- 
metry gives the other half of the parameter region. For 
small t/V and 0.15 < p < 0.85, the system exists in the 
solid state. From the above study on the system in the 
grand-canonical ensemble, this state can be regarded as 
a mixture of the solids with p = | and |. From p = 0.33 
to 0.67 and 0.1 < t/V < 0.2, the SS appears as in the 
grand-canonical ensemble. This result seems to be in 
sharp contrast with the phase diagram in the 2D trian- 
gular lattice at T = 0, in which the SS appear even for 
very small t/V^. This apparent discrepancy is closely 
related with how the SS evolves as T is raised. To see 
this, we study the system with t/V = c\/c2 = 0.1 and 
p = 0.5 by varying the value of /3. See FiglHl in which the 
specific heat C is shown as a function of C2. The result 




FIG. 8. Phase diagram in the canonical ensemble for C2 
50.0. 
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FIG. 9. Specific heat C for t/V = 0.1 and p = 0.5 as a 
function of C2. Tfiere exist two peaks indicating two second- 
order pfiase transitions. 
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lattice at T = 0^. In the grand-canonical ensemble, the 
most stable state appears for each value of the chemi- 
cal potential, and an inhomogeneous state is sometimes 
difficult to be observed. 

It is interesting how the PS state is realized in the 
present stacked triangular lattice^. In the 2D triangular 
lattice, the PS state is composed of domains of the solid 
and those of the SF. As shown in Fig(TOl the density and 
HC boson correlations for ci = 7.0 seem to imply that the 
SS is realized there, i.e., the solid-like behavior of n(r) 
and the SF-like behavior of G{r). Then we studied the 
correlations in the layered direction, and show the results 
in FigdU The results indicate that SF layers and solid 
layers appear separately as depicted in Fig(TT]in contrast 
to the 2D case. 
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FIG. 10. Number and HCB correlation functions in the 
canonical ensemble for C2 = 50.0 and p = 0.25. 

shows that there exist two phase transitions at C2 — 8.0 
and 22 both of which are of second order. By studying 
the correlation functions, we found that as T is raised, 
the phase transition from the SS to the solid state takes 
place first at C2 ^ 22 and then at C2 — 8.0 to the state 
without any long-range orders. See also Figs(T] and [5l 
Therefore the obtained finite-T phase diagram is consis- 
tent with that at T = in 2D. It should be remarked that 
in the 2D system at finite-T with p = 0.5, corresponding 
phase transitions take place almost simultaneously^. In 
Fig lTOl we show the correlation functions n(r) and G{r) 
for C2 = 50.0 and p = 0.25 in the canonical ensemble. 

In the phase diagram in the canonical ensemble shown 
in FiglHl the phase-separated (PS) state also exists, which 
does not appear in the grand-canonical ensemble. Similar 
result was obtained for the system of the 2D triangular 



FIG. 11. (Left) Number and HC boson correlation functions in 
the layered direction in the canonical ensemble for C2 = 50.0 
and p = 0.25. (Right)Sketch of a PS state in 3D. 

In conclusion, we have studied finite-T phase diagram 
of the HC boson system in the stacked triangular lattice. 
There exist various phases in the phase diagram and the 
obtained results are consistent with the phase structure 
of the HC boson system in the triangular lattice at T = 0. 
We also found that the SS first loses the superfluidity and 
then the density-wave order as T is raised. 
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